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ABSTRACT 

We identify the pressure-broadened red wings of the saturated potassium resonance lines at 7700 A as 
the source of anomalous absorption seen in the near-infrared spectra of Gliese 229B and, by extension, 
of methane dwarfs in general. In broad outline, this conclusion is supported by the recent work of Tsuji 
et al. 1999. The WFPC2 / band measurement of Glicsc 229B is also consistent with this hypothesis. 
Furthermore, a combination of the blue wings of this K I resonance doublet, the red wings of the Na D 
lines at 5890 A, and, perhaps, the Li I line at 6708 A can explain in a natural way the observed WFPC2 
R band flux of Gliese 229B. Hence, we conclude that the neutral alkali metals play a central role in the 
near-infrared and optical spectra of methane dwarfs and that their lines have the potential to provide 
crucial diagnostics of brown dwarf properties. 

The slope of the spectrum from 0.8 /imto 0.9 /imfor the Sloan methane dwarf, SDSS 1624+00, 
is shallower than that for Gliese 229B and its Cs lines are weaker. From this, we conclude that its 
atmosphere is tied to a lower core entropy or that its K and Cs abundances are smaller, with a preference 
for the former hypothesis. We speculate on the systematics of the near-infrared and optical spectra of 
methane dwarfs, for a given mass and composition, that stems from the progressive burial with decreasing 
T e ff of the alkali metal atoms to larger pressures and depths. 

Moreover, we surmise that those extrasolar giant planets (EGPs) that achieve T e fjsin the 800-1300 
K range due to stellar insolation will show signatures of the neutral alkali metals in their albedo and 
reflection spectra. We estimate that, due predominantly to absorption by Na D lines, the geometric 
albedo of the EGP r Boo b at A = 0.48 /mi is < 0.1, consistent with the new (and low) upper limit of 
0.3 recently obtained by Charbonneau et al. (1999). 

Subject headings: brown dwarfs, Gliese 229B, SDSS 1624+00, methane dwarfs, T dwarfs, extrasolar 
giant planets, alkali metals, atmospheres, non-gray spectral synthesis 



1. INTRODUCTION 

The discovery of Gliese 229B in 1995 was a milestone 
in the study of brown dwarfs, providing the first bona fide 
object with an effective temperature (T c g ~ 950 K) that 
was unambiguously substellar (Nakajima et al. 1995; Op- 
penhcimer et al. 1995). It also validated in dramatic fash- 
ion the allied theoretical and observational efforts of the 
many groups around the world that were engaged in brown 
dwarf research and started a race to discover and under- 
stand brown dwarfs that shows no sign of abating. Indeed, 
seven similar substellar mass objects, the so-called "T" 
or "methane" dwarfs, have since been discovered in the 
field by the Sloan Digital Sky Survey (Strauss et al. 1999; 
Tsvetanov et al. 1999), by the 2MASS survey (Burgasser 
et al. 1999), and by the NTT/VLT (Cuby et al. 1999). Su- 
perficially, these all seem to be clones of Gl 229B, but 
there are subtle spectral differences that no doubt reflect 
true differences in mass, metallicity, and effective tempera- 
ture. In addition to this flurry of brown dwarf discoveries, 
the first new stellar spectral class in almost 100 years, the 
"L type," has been defined, with ~100 members having 
been found to date by 2MASS (Kirkpatrick et al. 1999) and 
DENIS (Delfosse et al. 1997; Tinney et al. 1998). These 
objects are characterized by the clear onset of refractory 
heavy metal depletion, in particular that of titanium and 
vanadium, and by the growth in strength of the alkali 



metal lines of cesium, potassium, rubidium, and sodium 
in their optical and near-infrared spectra. L dwarfs also 
show distinct bands of FeH and CrH. Such transitions 
in chemical makeup with decreasing T e g are in keeping 
with theoretical predictions of the molecular consistuents 
of substellar atmospheres (Burrows and Sharp 1999; Feg- 
ley and Lodders 1996) that can serve to anchor and define 
the sequence of L dwarf spectral subtypes (Kirkpatrick 
et al. 1999). The L dwarfs have T off sfrom ~1500 K to 
^2000 K and constitute the spectroscopic link between M 
dwarfs and Gl 229B-like objects (Kirkpatrick et al. 1999). 
Many L dwarfs are also brown dwarfs (with masses below 
the stellar boundary at ^0.075 M©), but some will turn 
out to be stars very near the stellar edge. The current 
ambiguity reflects the newness of this subject. 

In this paper, we address T (methane) dwarf spectra in 
the optical and near infrared, focussing predominantly on 
Gl 229B as the prototype. Marley et al. (1996) and Allard 
et al. (1996) analyzed the full spectrum of Gl 229B (Op- 
penheimer et al. 1995; Geballe et al. 1996; Oppenheimer 
et al. 1998; Schultz et al. 1998) and obtained reasonable 
fits from the J band (~ 1.2 /jm) through the N band 
(~ 10 /mi). They concluded that its T cff was -900-1000 
K and that its gravity (g) was — 3 x 10 4 — 10 5 cm s~ 2 . 
The gravity error bars were large and translated into a 
factor of ~2-3 uncertainty in its inferred mass and age. 
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They also concluded that the atmosphere of Gl 229B is in- 
deed depleted in the refractory heavy elements, such as Al, 
Si, Ca, Fe, and Mg, just as the theory of condensate for- 
mation and rainout in cool molecular atmospheres would 
suggest, and that its spectrum longward of — 1 /imcan be 
fit, in the main, by a simple mixture of H2O, CH4, and 
H 2 . (Noll, Geballe, and Marley (1997) have detected CO 
at 4.67 /xmin Gl 229B as well.) However, neither Mar- 
ley et al. (1996) nor Allard et al. (1996) were able to fit the 
near-infrared observations between 0.8 /imand 1.0 /imand 
the theoretical excesses in flux ranged from 10 to 100. The 
line marked "Clear" in Figure 1 demonstrates the typical 
discrepancy between theory, that otherwise fits rather well 
at longer wavelengths, and the observed spectrum of Gl 
229B (Lcggett et al. 1999). Clearly, some ingredient is ab- 
sorbing in the blue, creating a very steep spectrum with a 
deep red (infrared) cast. 

Attempts to fit this problematic spectral interval were 
made by Golimowski et al. (1998), who invoked a Tsuji 
et al. (1996) spectrum that retained TiO in the atmosphere 
to take full advantage of TiO's strong absorption in the 
blue, and Griffith et al. (1998), who hypothesized that a 
population of small photochemical haze particles analo- 
gous to the red Titan Tholins (Khare and Sagan 1984) 
resided between -1400 K and -1800 K in Gl 229B's atmo- 
sphere. Indeed, TiO does give roughly the correct contin- 
uum slope, but it imposes the characteristic TiO bands on 
the spectrum that are not seen in either late L or T dwarfs. 
Furthermore, chemical abundance studies show that Ti 
and V are depleted (rained out) from the atmosphere to 
depths below the silicate clouds and to temperatures near 
2000 K, far below the Gl 229B's visible atmosphere. As to 
the suggestion by Griffith et al. , the dependence on wave- 
length between 0.85 /imand 0.93 /zmof the imaginary in- 
dex of refraction of the hazes inferred by construction was 
— 5x steeper than the corresponding published Tholin or 
polyacetylene values. Moreover, though Gl 229B has in 
its primary, Gliese 229A, a nearby source of UV photons 
needed to generate a haze, the late L dwarfs and the new 
field methane dwarfs do not, yet they demonstrate similar 
profound absorption in the blue. 

2. THE PIVOTAL ROLE OF THE K I RESONANCE 
DOUBLET AT ~7700 A 

What motivated Griffith et al. to posit the presence of 
a photochemical haze was the need for a continuum ab- 
sorber to suppress the "blue" flux and to act shortward of 
—1.0 /imdown to at least 0.8 /xm (Figure 1). Small Mie 
scattering particles might quite naturally have fit the bill. 
However, there is another, ready-made, source of contin- 
uum opacity, the red wings of the K I resonance lines at 
7665 A and 7699 A, that are in just the right place, short- 
ward of the Gl 229B data, with what appear to be just 
the right strength and opacity slope, to fully explain the 
anomalous data. The K I resonance doublet, due to the 
As 2 Si/2 — 4p 2 Pi/2,3/2 transitions, comes into its own only 
if other sources of opacity that dominate in M dwarfs, such 
as TiO and VO, are absent. The formation of refractories 
and their subsequent rainout accomplishes just that. 

The concept of "rainout" warrants further explanation: 
When condensates form, there are two possible general 
categories for their ultimate distribution within an atmo- 
sphere. The condensate (solid grain or liquid drop) can 



remain well-mixed with the atmosphere above the con- 
densation level, or the condensate can rain out of the at- 
mosphere. In the latter case, the exact vertical profile of 
the condensate depends upon poorly understood micro- 
physical processes. However, guided by experience with 
planetary atmospheres, we expect the condensate to form 
a cloud layer of some finite thickness. Both the condcnsiblc 
gas and its condensate are thus depleted above the cloud 
top. We refer to this process as "rainout." In the other, 
well-mixed, case the condensate can continue to react with 
the gas; we refer to this case as complete thermochemical 
equilibrium. Based on the Gliese 229B data, on experi- 
ence with the planets of our solar system, and on general 
physical grounds, some form of rainout seems to occur. 

As shown by Burrows and Sharp (1999), Fegley and 
Lodders (1996), and Lodders (1999), the alkali metals are 
less refractory than Ti, V, Ca, Al, Fe, and Mg and survive 
in abundance as neutral atoms in substellar atmospheres 
to temperatures of 1000 K to 1500 K. This is below the 
1500 K to 2500 K temperature range in which the sili- 
cates, iron, the titanates, corundum, and spinel, etc. con- 
dense. Hence, in the depleted atmospheres of cool brown 
dwarfs alkali metals quite naturally come into their own. 
Figures 2 and 3 show some representative abundance pro- 
files for the alkali metals Li, Cs, K, and Na, with and 
without rainout. A T cff =950 K and g=10 5 cm s~ 2 Gl 
229B atmosphere model from Burrows et al. (1997) and 
Burrows and Sharp (1999) was used. As demonstrated in 
Figures 2 and 3, with or without rainout elemental potas- 
sium and sodium persist to low temperatures to near the 
top of the T dwarf atmosphere. However, the distribu- 
tions of the atomic alkalis are not the same, with Li and 
Cs being deeper than Na and K. At lower temperatures 
in an atmosphere, the atomic forms give way to the chlo- 
rides and hydroxides {e.g., LiOH). With rainout (Figure 
2), below —1000 K both sodium and potassium exist as 
sulfides (Na 2 S and K 2 S) (Lodders 1999). Without rainout 
(Figure 3), a situation we deem unlikely, complete chemi- 
cal equilibrium at low temperatures requires that sodium 
and potassium reside in the feldspars, high albite and sani- 
dine. If such compounds formed and persisted at altitude, 
then the nascent alkali metals would be less visible, par- 
ticularly in T dwarfs. For either Figures 2 and 3 
demonstrate, all the elemental alkalis reside in the lower- 
pressure/lower-temperature reaches of cool, substellar at- 
mospheres. 

Lines of Cs I at 8521 A and 8943 A, Rb I at 7800 A 
and 7948 A, Li I at 6708 A, and Na I at 8183/8195 A have 
already been identified in L dwarf spectra. In addition, as 
Figure 4 and Kirkpatrick et al. (1999) demonstrate, the K 
I doublet is clearly dominant in the spectra of such late 
L dwarfs as 2MASS-1228, 2MASS-0850, Dcnis-0205, and 
2MASS-1632. Therefore, it is quite natural to explain the 
steep red slope of the new T dwarfs as being caused by the 
red wing of a saturated and pressure-broadened K I fea- 
ture. This is the thesis of our paper. Moreover, Gl 229B's 
I band flux (Mj —20.76), as measured using WFPC2 on 
HST (Golimowski et al. 1998) can be explained by the 
same K I resonance feature and the Gl 229B R band flux 
(Mr — 24.0), also measured by Golimowski et al. (1998), 
can be explained by a mix of the sodium D lines at 5890 
A, the lithium line at 6708 A, and the blue tail of the K 
I resonance doublet. Note that the 1.25 iim lines of ex- 
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cited K I have also been identified in T dwarfs (Strauss 
et al. 1999; Tsvetanov et al. 1999), so that the presence of 
potassium, at least at the higher temperatures in their at- 
mospheres, is in no doubt. The alkali metal lines may well 
prove to be the key to probing the atmospheric structure of 
the methane dwarfs, since their distributions and relative 
depths (c/., Figure 2) will be reflected in the systematic 
dependence of T and L dwarf spectra with T e s and g. 

3. PROFILES OF THE ALKALI METAL RESONANCE LINES 

Before we proceed with a discussion of our synthetic Gl 
229B spectra between 0.3 /xmand 1.5 /imwe discuss our al- 
gorithm for handling the pressure broadening of the alkali 
metal lines by molecular hydrogen. The line lists contain- 
ing the wavelength of the transition, the lower excitation 
energy, the log (5/) value, and the quantum numbers of the 
participating states were obtained from the Vienna Atomic 
Line Data Base (Piskunov et al. 1995). The major transi- 
tions of immediate relevance are those that correspond to 
the Na D lines at 5890 A and the K I resonance lines at 
7700 A. Given the high H 2 densities in brown dwarf at- 
mospheres, the natural widths (for Na D, ~0.12 mA) and 
Dopplcr widths of these lines are completely overwhelmed 
by collisional broadening. However, in general the line 
shapes are determined by the radial dependence of the 
difference of the perturber/atom potentials for the lower 
and upper atomic states (Griem 1964; Breene 1957,1981) 
and these are rarely known. The line cores are determined 
by distant encounters and are frequently handled by as- 
suming a van der Waals interaction potential with an adi- 
abatic impact theory (Wcisskopf 1933; Ch'en and Takeo 
1957; Dimitrijevic and Peach 1990) and the line wings are 
determined by close encounters and are frequently han- 
dled with a statistical theory (Holtzmark 1925; Holstcin 
1950). The transition between the two regimes is near 
the frequency shift (Ai/), or detuning, associated with the 
perturbation at the so-called Weisskopf radius (p w ), from 
which the collision cross section employed in the impact 
theory is derived (Spitzer 1940; Anderson 1950). In the 
simple impact theory, the line core is Lorentzian, with a 
half width determined by the effective collision frequency, 
itself the product of the perturbcr density, the average rel- 
ative velocity of the atom and the perturbcr (v), and the 
collision cross section (np^). If the frequency shift (Af) 
due to a single perturber is given by C n /r n , where r is the 
intcrparticle distance, then p w is determined from the con- 
dition that the adiabatic phase shift, J_ 2irAvdt 1 along 
a classical straight-line trajectory, with an impact param- 
eter p w , is of order unity. This yields p w oc (C n /v) 1 ^ n ^ 1 \ 
For a van der Waals force, n = 6. In the statistical the- 
ory, the line shape is a power law that goes like 1/Av~^~, 
and this is truncated (cut off) by an exponential Boltz- 
mann factor, e~ v °( Ts ^ kT , w j iere Vb(r s ) is the ground-state 
perturbation at the given detuning. The detuning at the 
transition between the impact and statistical regimes is 
proportional to {v 6 /C n )°- 2 (Holstein 1950). 

All this would be academic, were it not that for the 
Na/H2 pair the simple theory in the core and on the red 
wing is a good approximation (Ncfcdov, Sinel'shchikov, 
and Usachev 1999). We use this theory here. For the 
Na D lines perturbed by H 2 , we obtain from Nefedov, 
Sinel'shchikov, and Usachev (1999) a C 6 of 2.05 x 10~ 32 



cgs and a transition detuning, in inverse centimeters, of 30 
cm~ 1 (T/500.K') ' 6 , where T is the temperature. For the K 
I resonance lines, we scale from the Na D line data, using 
a Cq of 1.16 x 10~ 31 cgs, itself obtained from the theory of 
Unsold (1955). This procedure yields a transition detun- 
ing for the 7700 A doublet of 20 cm" 1 (17500^ ) ' 6 . From 
Nefedov, Sinel'shchikov, and Usachev (1999), we see that 
for a variety of perturbing gases the exponential cutoff for 
the Na D lines can be (for temperatures of 1000-2000 K) a 
few xlO 3 cm" 1 . The difference between 5890 A and 7700 
A, in inverse wavenumbers, is ~4000 cm -1 and that be- 
tween 7700 A and 1.0 /mi is only 3000 cm -1 . Hence, it is 
reasonable to expect that the dctunings at which the line 
profiles are cut off can be much larger than the Lorentzian 
widths or the impact/statistical transition detunings of 
tens of cm -1 . Since we as yet have no good formula for 
the exponential cutoff term, we assume that it is of the 
form e -i hAl/ /kT ^ wnere q j s an unknown parameter. Com- 
paring with the examples in Nefedov, Sinel'shchikov, and 
Usachev (1999), q may be of order 0.3 to 1.0 for the Na/H 2 
pair. Without further information or guidance, we assume 
that it is similar for the K/H 2 pair. We stress that this 
algorithm is merely an ansatz and that a more compre- 
hensive theory based on the true perturber potentials is 
sorely needed. Nevertheless, whatever the detailed line 
shape, as we will show, the basic conclusion that the K I 
resonance doublet is the "mystery" absorber in the near- 
infrared spectra of T dwarfs seems robust. 

Figure 5 depicts our opacity spectrum versus wavelength 
for the K I doublet at T = 2000 K and 1 bar pressure, for 
three different parametcrizations. Included are other, non- 
resonant, potassium lines excited by the high temperature. 
The solid line is for the standard Lorentzian theory, with- 
out a transition to the 3/2-power law. The other lines 
depict the line profile, corrected in the wings and with an 
exponential cutoff, using either q = 1.0 or q = 0.5. Note 
that for q = 1.0, beyond the core the transition to the 3/2- 
power law from the Lorentzian "2" -power law results in a 
larger opacity, which is then brought below the Lorentzian 
profile by the exponential. The same behavior is seen for 
q = 0.5, but the exponential brings the profile below the 
Lorentzian at a larger detuning. In all cases, the slope red- 
ward of the 7700 A doublet between 0.8 /mi and 1.0 /xmis 
of just the sort to affect the desired reddening. 

4. MODEL SPECTRA FOR GLIESE 229B 

Gliese 229B has been well-studied from the red to ~10 
pm and has a well-measured distance (5.8 parsecs) . Hence, 
with moderate resolution spectra in hand and absolute 
fluxes it is natural to use Gl 229B to demonstrate our 
basic thesis about the near-infrared and optical spectra 
of the entire T dwarf class. On Figure 1, in addition to 
depicting the "Clear" atmosphere fluxes from 0.5 pmto 
~1.45 pm , we give the synthetic spectra for four differ- 
ent models of Gl 229B that include the lines of the al- 
kali metal atoms, in particular the K I doublet at 7700 
A, employing the formalism of § 3. The Gl 229B data 
are taken from Leggett et al. (1999). Rainout abundance 
profiles of the alkali metals were used (c/., Figure 2). As 
Figure 1 demonstrates, the fit between 0.85 /mi and 1.0 
/mi is quite good and does not seem to require a layer of 
red particulates. We were guided in our choices of T e ffS , 
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gravities, and metallicities for these models by the work 
of Marley et al. (1996), Allard et al. (1996), and Griffith 
et al. (1998). Three of the models, those with higher fluxes 
around 0.7 /jmand lower fluxes near 1.0 /xmused a q of 0.4 
and assumed alkali metal abundances of 0.3 times solar 
(Anders and Grevesse 1989). They had T e g/g pairs of 
[900 K/10 5 (cgs)], [950 K/10 5 (cgs)}, and [800 K/3 x 10 4 
(cgs)]. From Burrows et al. (1997), the ages and masses for 
these models are 1.46 Gyr/34.9 Mj , 1.26 Gyr/35.3 Mj , 
and 0.44 Gyr/14.8 Mj , respectively. The fourth model, 
that with the lowest fluxes near 0.7 /xmand the highest 
fluxes near 1.0 /jm , employed a q of 1.0, had a T e gf /g pair 
of [950 K/10 5 {cgs)], and assumed that the alkali metal 
abundances were solar. Note that values for q of 0.4 and 
1.0 are merely possibilities and that wc do not necessar- 
ily advocate either one. As stressed in §3, the far wings 
of the neutral alkali metal lines, perturbed by molecular 
hydrogen, have not been properly calculated. In fact, to 
our knowledge there is no other environment in which the 
neutral alkali line strengths at 1000-3000 A detunings have 
ever before been needed. 

Figure 6 depicts the behavior of the theoretical spec- 
tra as a function of q, from to 1, for the [950 K/10 5 
(cgs)] model (with 0.3xsolar alkali metal abundances). 
Also included is the predicted spectrum for the uncor- 
rected Lorentzian K I line profile (dotted). Though Figure 
6 depicts a range of near-infrared and optical spectra, the 
character of these spectra is similar. However, as a conse- 
quence of the current ambiguity in the proper shape of the 
7700 A profile, which we have parameterized with q, we 
can not yet tie down the alkali metal abundances to better 
than perhaps a factor of three. Nevertheless, the good fits 
indicated in Figure 1 are quite compelling. Rather stun- 
ningly, both the measured WFPC2 I band flux near 0.814 
/xmand the WFPC2 R flux near 0.675 ^m (Golimowski 
et al. 1998) fit as well. Hence, we conclude that the I 
band flux is determined by the red wing of the K I reso- 
nance doublet and that the R band is a consequence of the 
red wings of the Na D lines (~ 5890 A), the blue wings 
of the K I doublet, and, perhaps, the lithium line at 6708 
A. The broad width of the / band filter extends it to be- 
yond the saturated center of the doublet at 7700 A. Our 
theoretical models have absolute / band magnitudes from 
21.0 to 21.3, compared with the measured value of ^20.76. 
Given the calibration problems when an underlying spec- 
tral slope is so steep and so unlike that of the standard 
stars, we conclude that the fit is indeed good. As Figure 1 
demonstrates, our R band predictions bracket the R band 
measurement. 

We emphasize that one should be very careful about 
interpreting band magnitudes as fluxes. There are two 
points to make in this regard: 1) bands are spread over 
a sometimes broad range of wavelengths and, 2) magni- 
tudes are calibrated on stars with certain underlying spec- 
tra that may not be similar to the spectrum of the object 
under study. If these spectra are very different, large er- 
rors can be made in assigning magnitudes, as well as in 
assigning fluxes at a given wavelength in the band. A 
quick look at the spectrum of Gl 229B from 0.8 to 1.0 /xm 
will show why this should be of concern; the calibration 
stars can not have such steep slopes and the errors in the 
measured magnitudes might be very large, many tenths 
of magnitudes. This, and extravagant band width, are 



particularly relevant for the Gunn r and i band measure- 
ments of Matthews et al. (1996). Given this, we compare 
with just the WFPC2 R and I band data of Golimowski 
et al. as a set from a single group and depict the breadth 
of these bands on our figures. 

We note in passing that the discrepancy between the 
depth of the predicted and measured flux trough between 
the Z (~1.05 yum) and the J (~1.25 |«m) bands may be 
resolved with a lower abundance of water and, hence, oxy- 
gen (Griffith et al. 1998) and that the shape of the J band 
can not be fit without the methane bands at ~1. 15 — 1.2 
/imand ~ 1.3 /xmto sculpt it. Hence, the shape of the J 
band in T dwarfs is not just a consequence of bracketing by 
H 2 absorption features, but requires methane and is an- 
other signature of methane is their atmospheres. However, 
this point, as well as the inferred elemental abundances, 
are not the subject of the current paper and we defer a 
fuller discussion to a later work. 

The variety of models that fit the Gl 229B data, given 
the current state of the theory and observations, was noted 
in Marley et al. (1996) and Allard et al. (1996) and re- 
flects the fact that different models with different T e ff / g 
pairs can have similar atmospheric temperature-pressure 
profiles and luminosities. The atmospheres of substcllar 
objects are convective at depth and radiative on the pe- 
riphery. For a given composition, brown dwarfs and ex- 
trasolar giant planets (EGPs; Burrows et al. 1995) are a 
two-parameter family; given two independent quantities 
such as T e ff , luminosity, radius, gravity, or age, one can 
derive, with theory, any of the others. In particular, as dis- 
cussed in Hubbard (1977) and Saumon et al. (1996), in an 
approximate sense, the entropy (S) at depth is a function 
of T e ff and g and this entropy determines the tempera- 
ture/pressure profile. (This is not to say that the atmo- 
spheres are adiabatic, merely that the T/P profile of the 
outer radiative zone can be calculated from S, T e g , and g.) 
Hubbard (1977) showed that S is approximately a function 
of the combination T^/g 1 / 6 . From Marley et al. (1996), 
we obtain rough power-law relations for mass (M) , radius 
(R), and age (t), as a function of g and T c ff : 

/ a \ 0.64 /7 1 x 0.23 

» = (w) 

/ 10 5 \ 0-18 /T „\ 0.11 

»- 6.7x10^) @) 

t-i-Mw) (55) • (1) 

where g is in cm s~ 2 . From these equations, and the depen- 
dence of S on T ff and g cited above, we derive approxi- 
mately that, for a given core entropy and, hence, for a 
given T/P profile, M ~ t - 55 . Using the same set of power 
laws, we can derive that, for a given luminosity, M ~ i - 43 . 
Hence, since these two power-law relations are similar, fits 
to Gl 229B roughly define a trajectory in M-t space with 
a ~ 0.5 power. This is why low- mass/short- age models 
and high-mass/long- age models both fit Gl 229B, at least 
for the low-resolution comparisons thus far attempted. In 
principle, this degeneracy can be broken once the the- 
ory and the observations are better constrained. Figure 7 
shows temperature/pressure profiles for a few models from 
Burrows et al. (1997). The higher curves are for higher 
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T c ffsand lower gravities. Note that the lower-T e ff , lower- 
g model [800 K/3 x 10 4 (cgs)] is in the same vicinity as 
the higher-Tcff , higher-^ model [900 K/10 5 (cgs)}, both of 
which are good fits to Gl 229 B (Figure 1). Note also that 
the [1100 K/5 x 10 4 (cgs)}, [400 K/5 x 10 4 (cgs)}, and [800 
K/10 5 (cgs)] models depicted for comparison do not fit Gl 
229B, as we might have surmised from Figure 7. 

5. THE SLOAN DWARF: SDSS 1624+00 

The models and data depicted in Figure 1 capture the 
Cs lines at 8521 A and 8943 A, indicating that cesium is 
in Gl 229 B's atmosphere in modest abundance between 
1200 K and 1400 K. The depths to which cesium can be 
found are shown in detail for a specific model in Figures 2 
and 3 and approximately for the variety of models by the 
position of the intercept in Figure 7 of the "Cs/CsCl=l" 
line with a given T/P profile. Similarly, the depths of 
atomic K, Na, Rb, and Li can be found. Shown on Fig- 
ure 7 along with the Cs/CsCl=l line is a corresponding 
complete equilibrium (no rainout) "K/KC1=1" trajectory. 
As Figure 2 demonstrates, with rainout (a more realistic 
situation), atomic potassium persists to much lower tem- 
peratures, perhaps shifting the K/KC1=1 line on Figure 
7 down by 200-300 K. Nevertheless, qualitatively and im- 
portantly, the lower the T/P profile in Figure 7, the deeper 
are the atomic cesium and potassium. Depth and higher 
pressure imply a higher column density of H 2 0, CH4, and 
H2 and, hence, a higher optical depth of these competing 
absorbers. Therefore, for a given composition, at lower 
T e fi and higher g, both potassium and cesium arc more 
deeply buried than at higher T e ff and lower g. The ex- 
pected upshot is the gradual diminution with decreasing 
core entropy of the strengths of the cesium and potassium 
features. A decrease in the effect of the K I resonance dou- 
blet at 7700 A would result in an increase in the influence 
below 1.0 /imof the H2O, CH4, and H2 features that are 
not as steep as the red wing of the doublet. As a result, the 
spectrum short of ^0.95 /im would shallow and the 7700 A 
region of the spectrum would fill in. The Rb lines at 7800 
A and 7948 A may make an appearance in the K I trough, 
but they will still be rather weak. Note that, for a given 
T e g / g pair, lowering the abundances of the alkali metals 
would raise the flux in the 0.7-0.95 /im region, but it might 
also make the Z — J color bluer than observed for the cur- 
rent crop of T dwarfs. Furthermore, the presence of the 
potassium line at 1.25 /imboth in Gl 229B and in the two 
Sloan dwarfs makes this explanation suboptimal. Never- 
theless, potential variations in abundances and metallicity 
from object to object do complicate the analysis. 

The shallowing of the spectrum shortward of 0.9 /im and 
the weakening of the cesium features at 8521 A and 8943 
A are indeed seen in the preliminary spectrum of SDSS 
1624+00 (Strauss et al. 1999). Strauss et al. note that the 
spectra of Gl 229B and SDSS 1624+00 seem uncannily 
similar and conclude that they must have similar mass, 
age, and T e ff . Given the lower apparent flux of SDSS 
1624+00, they put this object further away than Gl 229B, 
at ~10 parsecs. However, theory implies the J — H and 
J — K colors don't change hugely in the 600 K to 1200 
K range (Burrows et al. 1997). As a consequence, for a 
given gravity or mass, a decrease in T e g merely, though 
approximately, translates the absolute spectrum down to 
lower fluxes. Hence, there is a need for more subtle di- 



agnostics of a substellar object's properties. Such diag- 
nostics may be the strengths of the alkali metal lines in 
the optical and near infrared. Figure 8 compares the ob- 
served SDSS 1624+00 spectrum with a theoretical model 
spectrum, as well as with the observed Gl 229B spectrum. 
The model (dotted) has an atmosphere with T c ff =700 K 
and 5=3 x 10 4 cm s~ 2 . Its core entropy is below those 
of the Gl 229B models. We theorize that the shallower 
near-infrared spectrum of SDSS 1624+00 with respect to 
Gl 229B, as well as the weakness of the Cs lines, are both 
consequences of a lower core entropy T/P profile. This 
could mean a lower T c ff , a higher g or mass, or some suit- 
able combination that buries the tops of the potassium 
and cesium distributions at higher pressure and, hence, 
larger H2O optical depths than found in Gl 229B. How- 
ever, given the remaining theoretical ambiguities, we can't 
completely rule out the possibility that lower K and Cs 
abundances are the explanation. Note that it is still pos- 
sible that T off for SDSS 1624+00 is much larger than 700 
K. However, given its provisional cesium line strengths, we 
would then require its gravity, and presumably its mass, 
to be correspondingly higher to achieve the lower core en- 
tropy we suggest can explain its spectrum. Importantly, if 
improved spectra were to show that Cs is much stronger in 
SDSS 1624+00 than its preliminary, low-resolution spec- 
trum now indicates, a higher T e g , perhaps higher than 
that of Gl 229B, would be indicated. This emphasizes the 
importance of the alkali metal lines as diagnostics of the 
T dwarfs and the need to obtain better spectra. Note that 
with the [700 K/3 x 10 4 cgs] model, the absolute Z and J 
band fluxes of the SDSS 1624+00 are below those for Gl 
229B and match only near A = 0.85 /im. If the T c g of 
SDSS 1624+00 is indeed lower than that of G1229B, SDSS 
1624+00 could be closer than the 5.8 parsecs of Gl 229B. A 
smaller distance would increase the inferred number den- 
sity of SDSS 1624+00-like objects, perhaps beyond what 
is credible. Hence, obtaining a parallax for this new T 
dwarf will be a crucial prerequisite for future substantial 
improvements in characterizing it. 

Importantly, any T c ff estimate for SDSS 1624+00 must 
be confirmed by a comparison of models and observations 
at longer wavelengths, similar to that performed for Gl 
229B (Marley et al. 1996; Allard et al. 1996). As with 
Glicsc 229B, simultaneously achieving a satisfactory fit to 
the optical flux, the near-infrared windows, and absorption 
features will be challenging. In particular, if SDS 1624+00 
has a lower T e ff than Gliese 229B, the water and methane 
band depths will likely be deeper. Yet the methane opacity 
database at 1.7 /imis too poor to permit detailed compar- 
isons and the 2.3 /^mband is also influenced by H 2 opacity, 
the relative strength of which depends upon the metallic- 
ity. L' band spectra, which probe the methane fundamen- 
tal band at 3.3 ^m where there are no competing absorp- 
tions and the opacity database is better understood, may 
provide the best point of comparison. 

However, for a given mass and composition, the "his- 
tory" of the alkali metal lines in L and T dwarfs may 
follow the following sequence: At higher T c ff s (in the L 
dwarf range), the depiction of the refractory elements al- 
lows the strengths of the alkali metal lines to wax. The ce- 
sium lines peak before reaching the T e ff of Gl 229B, though 
they are still prevalent in such atmospheres, and potas- 
sium determines the character of the spectrum from 0.7 
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/xm through 1.0 fim . As T e g decreases further, the alkali 
metals are buried progressively more deeply and begin to 
wane in strength. Due to this burial of potassium, the 
slope of the spectrum from 0.8 /umto 1.0 fim decreases. 
Since the abundance distribution of the atomic form of 
each alkali metal is different (Figures 2 and 3), the T e ffsat 
which the lines of a given alkali start to wane will be dif- 
ferent, with the effects of Li, Cs, Rb, K, and Na waning in 
approximately that order with decreasing T e g (at a given 
mass and metallicity). The specific numbers will depend 
upon the actual effects of rainout on the alkali metal pro- 
files, but the systematics should not. Eventually, perhaps 
below T e ffS of 500-600 K, the effects of the neutral alkali 
metals are eclipsed, with clouds of H 2 eventually effect- 
ing the next important change in the character of brown 
dwarf atmospheres. Given this systematics, whatever the 
true T cff of SDSS 1624+00, we would predict that both 
above Gl 229 B's T e ff , near the T e ff juncture between the 
L dwarfs and the T dwarfs, and below a T e ff of ~800 K, 
the troughs around the K I resonance doublet at 7700 A 
would be partially filled and the slope of the "optical" 
spectra between 0.8 and 1.0 A would shallow. 

6. SUMMARY 

We conclude that the anomalous absorption seen in the 
near-infrared and optical spectra of all the T dwarfs dis- 
covered to date is due to the red wings of the saturated K 
I resonance lines at 7700 A. This theory also explains the 
WFPC2 I and R flux measurements made of Gl 229B, with 
the Na D lines at 5890 A helping to determine the strength 
of the R band. There are still ambiguities in the T c ff s , 
gravities, and compositions of Gl 229 B, in particular, and 
of T dwarfs, in general, but a sequence from Gl 229B to 
SDSS 1624+00 of decreasing core entropy is seen to be 
consistent with the expected systematics of the tempera- 
ture/pressure and alkali-metal abundance profiles. (How- 
ever, lower K and Cs abundances for the Sloan dwarf or 
a reassessment of its cesium line strengths once higher- 
resolution spectra become available can not yet be ruled 
out.) Silicate grains are expected at depths of 1500 K to 
2500 K. Their formation is inferred in the late M dwarfs 
(Jones and Tsuji 1997) and their eventual burial is seen in 
the appearance at lower T e ffSof the new L spectral class 
(Kirkpatrick et al. 1999). Complete refractory element de- 
pletion ushers in a phase of neutral alkali metal dominance 
at short wavelengths, which persists until the tops of the 
alkali metal distributions are buried at higher pressures 
and higher H 2 0, CH4, and H 2 optical depths. This oc- 
curs in low-entropy atmospheres, though low alkali metal 
abundances can mimic the same effect somewhat. We note 
that if the Lorentzian theory for the line shapes of the Na 
D lines were used, due to the high relative abundance of 
sodium, their influence would stretch 7000 cm -1 to 1.0 
/im and would flatten the top of the Z band, contrary to 
observation. This is indirect evidence for the action of an 
"exponential cutoff' on the Na D line shape. Even with 
such a cutoff, the K I resonance line dominates the spectra 
of Gl 229 B and its ilk from 0.7 /mi to 1.0 /jm . For Gl 229 
B, we predict the presence of a saturated absorption line 
around 7700 A, with higher fluxes on either side. We also 
expect a deep, broad trough around the Na D lines at 5890 
A (see Figure 1). 



At the time of this writing, we became aware of a paper 
by Tsuji et al. (1999) that has also concluded that the K 
I feature is the major, "mystery" absorber at short wave- 
lengths. We wholeheartedly concur with this conclusion. 
However, contrary to Tsuji et al. , we determine that dust 
may not be required to achieve a good fit shortward of 
1.0 /im. Moreover, we believe that the Gliese 229B / 
band flux is well-fit by the red wing of the K I doublet 
alone. In addition, if the Na D lines are included in a nat- 
ural way, an acceptable fit to the WFPC2 R band flux is 
achieved as a byproduct. However, filling in the troughs 
between the Z and J bands and between the H and K 
bands may still require some combination of grains, as 
Tsuji et al. suggest, and sub-solar metallicity, though sub- 
solar metallicity and improvements in the molecular opac- 
ity databases alone may be sufficient. If silicate dust is 
important, it could be because of upwclling from below, 
a very non-equilibrium process that is not easily modeled. 
In addition, optical constants, particle shapes, and parti- 
cle sizes would need to be known. There is no credible 
guidance concerning the particle size; the guessed average 
particle radius could off by one or two orders of magnitude. 
This ambiguity translates into a correspondingly large am- 
biguity in the dust opacity. As Griffith et al. (1998) demon- 
strate, a sufficiently arbitrary dust model can be found to 
fit the Gl 229B spectrum. Nevertheless, a possible model 
might be one in which a small population of dust from mi- 
nor species at low pressures warms the atmosphere near 1 
bar, thus decreasing the depths of the water and methane 
absorption bands which originate near this pressure, yet 
providing little opacity elsewhere. However, we reiterate 
that, given the new-found importance of the K I feature, 
there is much less reason to evoke a population of dust or 
grains to fit methane dwarf spectra. 

Note that in this paper we have considered only ob- 
jects in isolation. However, we expect that those EGPs 
that achieve T e ffsin the 800-1300 K range due to stellar 
insolation will also show signatures of the neutral alkali 
metals. Recently, Charbonneau et al. (1999) have put an 
upper limit of 0.3 to the geometric albedo at A = 0.48 
/imof the planet orbiting r Boo. This is below some pub- 
lished predictions (Marley et al. 1999), though not others 
(Seager and Sasselov 1998), and may indicate the presence 
in its atmosphere of sodium and absorption by the Na D 
lines (Sudarsky, Burrows, and Pinto 1999). Our prelimi- 
nary estimate of the geometric albedo at A = 0.48 /imof 
such a planet, due to the influence of the neutral alkali 
metals, is < 0.1. Since stellar insolation is bound to create 
hazes (Marley 1998), such as absorb in the blue and UV 
in Jupiter, detailed reflection spectra of r Boo and similar 
"close-in" EGPs will be needed to disentangle the relative 
contribution of the alkali metals to EGP albedos. 
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Fig. 1. — The log of the absolute flux (F„) in millijanskys versus wavelength (A) in microns from 0.5 /imto 1.45 fimior Gliese 229 B, 
according to Leggett et al. (1999) (heavy solid), and for four theoretical models (light solid) described in the text. Also included is a model, 
denoted "Clear" (dotted), without alkali metals and without any ad hoc absorber due to grains or haze. The horizontal bars near 0.7 /mi and 
0.8 /im denote the WFPC2 R and / band measurements of Golimowski et al. (1998). 
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Fig. 2. — The fractional abundances of different chemical species involving the alkali elements Li, Na, K and Cs for a Gliese 229B model, 
with rainout as described in Burrows and Sharp (1999). The temperature/pressure profile for a T e g =950 K and g = 10 5 cm s~ 2 model, 
taken from Burrows et al. (1997), was used. Each curve shows the fraction of the alkali element in the indicated form out of all species 
containing that element. All species are in the gas phase except for the condensates, which are in braces { and }. The solid curves indicate the 
monatomic gaseous species Li, Na, K and Cs, the dashed curves indicate the chlorides, the dot-dashed curves indicate the hydrides and the 
triple dot-dashed curve indicates LiOH. Due to rainout, at lower temperatures there is a dramatic difference with the no-rainout, complete 
equilibrium calculation (Figure 3); high albitc and sanidinc do not appear, but instead at a much lower temperature the condensate Na2S 
(disodium monosulfide) forms, as indicated by the solid line in the lower left of the figure. The potassium equivalent, K2S, also forms, but it 
does so below 1000 K and is not indicated here. The difference between this Figure and Figure 3 is that almost all the silicon and aluminum 
have been rained out at higher temperatures, so that no high albite and sanidine form at lower temperatures. 
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Fig. 3. — The fractional abundances of different chemical species involving the alkali elements Li, Na, K and Cs for a Gliese 229B model, 
assuming complete (true) chemical equlibrium and no rainout (disfavored). The temperature/pressure profile for a T e g =950 K and g = 10 
cm s~ 2 model, taken from Burrows et al. (1997), was used. Each curve shows the fraction of the alkali element in the indicated form out of 
all species containing that element, e.g., in the case of sodium, the curves labeled as Na, NaCl, NaH and NaAlSiaOg are the fractions of that 
element in the form of the monatomic gas and three of its compounds. All species are in the gas phase except for the condensates, which 
are in braces { and }. The solid curves indicate the monatomic gaseous species Li, Na, K and Cs and the two condensates NaAlSiaOg and 
KAlSi3C>8, i.e., high albite and sanidine, respectively, the dashed curves indicate the chlorides, the dot-dashed curves indicate the hydrides 
and the triple dot-dashed curve indicates LiOH. 
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Fig. 4. — The log of the absolute flux (F„) in millijanskys versus wavelength (A) in A from 6500 A to 10000 A for Gliese 229 B and 
representative L dwarfs from Kirkpatrick et al. (1999). Included are 2MASP-J0345 (L0), GD 165b (L4), 2MASSs-J0850 (L6), and 2MASSW- 
J1632 (L8), in order of decreasing general flux level. The L dwarf spectra were put on an absolute scale using the distances and apparent J 
band magnitudes in Kirkpatrick et al. and should be used with extreme caution. 
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Fig. 5. — The log of the photon-potassium cross sections per particle, in cm , versus wavelength (A) in (im. The focus is on the region 
around the K I resonance doublet at 7700 A and for specificity a temperature and pressure of 2000 K and 1 bar have been assumed. The 
solid line is the Lorcntzian theory without corrections in the broad wings at large detunings. The short-dashed and long-dashed lines are for 
the corrected theory, with qs of 0.5 and 1.0, respectively. The crucial fact is the useful slope between 0.8 /imand 1.0 fim. 
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Fig. 6.— Similar to Figure 1, using the T eff =950 K and g = 10 5 cm s~ 2 Gliese 229B model, but for qs of 0.0, 0.25, 0.5, and 1.0. The 
lower the q, the lower the curve. The dotted line is the theoretical spectrum using the Lorentzian profile. Alkali metal abundances of 0.3x 
solar were assumed. The heavy solid line is the measured spectrum of Gliese 229B and the solid horizontal lines are the WFPC2 Gliese 229B 
R and / fluxes. 
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Fig. 7. — The pressure (in Atm) versus the temperature (in Kelvin) for various brown dwarf atmospheres, as a function of T e ff and gravity 
and at solar metallicity. The three central solid lines are for 900 K, 950 K, and 1000 K models with a gravity of 10 5 cm s~ 2 . The two outliers 
are models at temperatures (1100 K and 400 K) and entropies that are too high and too low, respectively, to fit the Gliese 229 B spectral 
data. The long-dash/short-dash curves are models with T c ff = 800 K, the higher one at 3 X 10 4 cm s~ 2 fits, while the lower one at 10 5 cm s~ 2 
does not. Also included are the lines that demark the true-equilibrium, no-rainout temperature/pressure trajectories for which atomic K and 
gaseous KC1 and atomic Cs and gaseous CsCl have equal abundances. While the Cs/CsCl line will be little changed by rainout, the K/KC1 
line will be shifted to lower temperatures by perhaps 200-300 K. These lines serve to indicate the character of the burial of the elemental 
alkali metals with decreasing core entropy. As the core entropy decreases, the pressures and column depths at the intersection points steadily 
increase. 
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Fig. 8. — The log of the absolute flux (F„) in millijanskys versus wavelength (A) in microns for the first Sloan T dwarf, SDSS 1624+00 
(heavy solid line), from Strauss et al. (1999). The light solid line is Gliese 229B, while the dotted line is for a model with T e £f =700 K and 
g = 3 X 10 4 cm s — 2 . The core entropy of this theoretical model is lower than that for those that fit Gliese 229B. Both the model and the 
Sloan data have a shallower slope than Gliese 229B from 0.85 fimto 1.0 ^mand weak cesium lines. 



